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Abstract 



In this Letter, we find suitable potentials in the multiple scalar fields scenario by using 
the Noether symmetry approach. We discussed three models with multiple scalar fields: 
N-quintessence with positive kinetic terms, N-phantom with negative kinetic terms and 
N-quintom with both positive and negative kinetic terms. In the N-quintessence case, the 
exponential potential which could be derived from several theoretic models is obtained from 
the Noether conditions. In the N-phantom case, the potential ^(1 — cos (^/^^fy))) which 
could be derived from the Pseudo Nambu-Goldstone boson model, is chosen as the Noether 
conditions required. In the N-quintom case, we derive a relation DV^ q = —DV^ p between 
the potential forms for the quintessence-like fields and the phantom-like fields by using the 
Noether symmetry. 
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1 Introduction 



Scalar field theory which is related to particle physics has become the generic playground for 
building cosmological models, both in the early and late accelerating periods of our universe 
[U E] . Although the dynamics of these accelerations is likely to contain several scalar fields, 
it is normally assumed that only one of these fields remained dynamically significant for 
a long time. However, realistic theoretical models, embedded in grand unified or super 
symmetric theories, must necessarily be theories of multiple fields. The simplest multiple 
scalar fields scenario which we will consider is first originated from the assisted inflation 
scenario [3J. The essential point of this scenario is that inflation is not driven by any 
single field, but a collection of N fields. These fields have the same initial conditions and 
potentials. This idea can be applied in vector field models as well pi]. 

Meanwhile, the observations suggest the equation of state (EoS) parameter of dark 
energy is in the range of— 1.21 < u < —0.89 [5]. Since the quintessence type of matter 
could not give the possibility that oj < —1, the extended paradigms (e.g. phantom and 
quintom) are proposed. Phantom type of matter with negative kinetic energy has well- 
known problems, but, nevertheless, was implicitly suggested in cosmological models and 
have also been widely studied as dark energy. It is phenomenologically significant and 
worthy of putting other theoretical difficulties aside temporally. Then, it is natural to ask 
why don't we discuss the multiple scalar fields with different kinetic terms. According to 
the classification of the scalar fields^, we can discuss three types of fields in the simplest 
multiple scalar fields scenario, which are the quintessence type of fields with positive kinetic 
terms, the phantom type of fields with negative kinetic terms, the quintom type of fields 
with both positive and negative kinetic terms. In this Letter, we call them N-quintessence, 
N-phantom, N-quintom for convenience. 

But, as in the single scalar field case, we have to ask how to choose the potentials from 
the various models for those multiple scalar fields. In this Letter, we will deal with this 
problem of choice from a point of view of symmetry. The Noether symmetry has been 
revealed as a useful tool for finding out exact solutions in cosmology. This is an interesting 
method to select models motivated at a fundamental level. 

This Letter is organized as follows. In section El we introduce the multiple scalar fields 

4 The quintessence with positive kinetic term was proposed in Ref. [6] ; the phantom with negative kinetic 
term was suggested in Ref. [7] ; and the quintom with both positive and negative kinetic terms was proposed 
in Ref. [8]. 
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models. In section [HI the Noether symmetry approach will be introduced and applied to 
both N-quintessence and N-phantom cases to get exact solutions. In section HI we discuss 
the application of Noether symmetry approach to N-quintom case in connection with its 
solution. In section [5j we give out the evolution of our universe in N-quintessence and 
N-phantom cases. Finally, a short summary will be presented in section El 

2 N-quintessence, N-phantom and N-quintom Scalar 
Field Model Scenario 

As stated in the introduction, usually, only one scalar field is enough to accelerate the 
universe, but a single field is not natural. The application of the multiple scalar fields in 
cosmology should be seriously considered. Here, we assume that the geometry of space- 
time is described by the flat FRW (Friedmann-Robertson- Walker) metric which seems to 
be consistent with today's cosmological observations 

3 

ds 2 = -dt 2 + a 2 (t)J2( dxi ) 2 , (!) 

where a is the scale factor. After setting the number of the scalar fields as N, the action of 
the multiple scalar fields can be written as 



5*0 = / d 4 xy/^g 



16nG 

i=i 



(2) 



where e = 1 denotes the quintessence fields with the positive kinetic terms, e = — 1 denotes 
the phantom fields with the negative kinetic term. Meanwhile, as we consider both the 
vector fields and the matter in the system, the total action is 

Stot — Sep + S m , (3) 

where S m is the action for matter. The density of the matter can be expressed as p m = 
Pmo( a o/°) 37 ) where p m0 is an initial constant and < 7 < 2. Here, we limit our analysis to 
7 = 1 which corresponds to the pressureless matter with P m = 0. 

We assume the vector fields are non-interacting, their influences on each other are 
through their effects on the expansion. Considering all the scalar fields have the same 
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potentials and initial conditions, action ([2j) could be simplified as 



S, 



01 



d 4 Xy/—g 



R 



16ttG 



+ N\ e^--V(4>) 



(4) 



When e = 1, we call the related scenario N-quintessence. While e = — 1, we call the related 
scenario N-phantom. 

For the N-quintom case, we assume the fields with same kinetic terms have the same 
potentials and initial conditions, the action can be written as 

R 



S, 



62 



d XyJ—g 



16nG 



+ N a 



-V(<j> q ) +N P 



- V{(j> p ) 



(5) 



where <p q is the scalar field with the positive kinetic terms, N q is the number of the corre- 
sponding quintessence type fields; P is the scalar field with the negative kinetic terms, N p 
is the number of the corresponding phantom type fields. This paradigm has been proved 
of crossing uj^ = — 1 when N g = N p = 1 [Hj. 



3 The Noether Symmetry Approach in N-quintessence 
and N-phantom 

In the case of N-quintessence and N-phantom, we take the scale factor a and the scalar field 
as independent dynamical variables in the system which the action (jlj) represents. Then 
the configuration space could be chosen as Q = (a, 0), while the related tangent space is 
TQ = (a, 0,a,0). To study the symmetries of the space under consideration, we need an 
effective point-like Lagrangian for the model whose variation with respect to its dynamical 
variables yields the correct equations of motion. However, based on action (J4j) , it is proper 
to make the point-like Lagrangian as 

A = + C m = 3aa 2 - 4- (e^f ~ a 3 V(4>)) + (6) 

m P i \ z J m pi 

where the Planck mass is = (87rG) -1 , and the term corresponds to the effects from 

m P i 

matter. 

Therefore, the total energy of the system Ec^ , could be written in this way 

E Cl = - & = a 3 ( e -^f- + NV{(j>) + Pm0 a- 3 - 3m 2 pl H 2 ) . (7) 



If the above equation being considered as a constraint, with the vanishing of the "energy 
function" , it is just the Friedmann equation 



H 2 = 



eN(j> 2 



NV((j))+p m0 a 



(8) 



Furthermore, for a dynamical system, the Euler-Lagrangian equation is 

Based on the Lagrangian, in the N-quintessence and N-phantom case, the variable qi is a 
and 0, respectively. When q { = a, the Raychaudhuri equation could be gotten 

U P , ^ eNm P l i2 m ll 

H= ^-{P<P + + Pm) = ^0 ~ 2 Pm ' ^ 

where the energy density and the pressure of scalar fields are 

P0 = ^0 2 + An/(0), (ii) 
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eN 



P,= ^—<f> 2 -NV (0). (12) 
What is more, the equation of state could also be obtained 

a = «f/2-yW (13) 

e0 2 /2 + \/(0) 

Obviously, in the N-quintessence case where e = 1, ui^ > — 1; in the N-phantom case where 
e = —1, u;^ < —1. Both of them could not cross oj^ = —1, that is why we also consider 
N-quintom. In the case of = 0, the Euler-Lagrangian equation is the equation of motion 

+ 3H + eVl = 0, (14) 



where the prime means = dV/dcj). For the different value of e, the quintessence makes 
the fields roll down the potential, while the phantom makes them roll up. 

The above equations coincide with the results calculated from the Einstein equations, 
and prove that the point-like Lagrangian is consistent with the dynamical system. 
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As is well known in (TTJ [T2J, [13] , Noether symmetry approach is a powerful tool in finding 
the solution for a given Lagrangian. From this method, it is possible to obtain a reduction, 
and possibly get a full integration of the system, whenever the cyclic variable of the system 
is found. The key point related to the Noether symmetry is a Lie algebra presented in the 
tangent space. Following [TTJ [T2], [13] , for the Lagrangian ([6]), firstly we define the Noether 
symmetry induced by a vector X on the tangent space TQ = (a, 4>, a, 0) which is 

d d d-d 

where a and /3 are generic functions of a and </>. The Lagrangian is invariant under the 
transformation X if 

L C ^ 1 + < ^ a ^ 1 _|_ -|- dfidd Q Qg\ 

da dt da da dt da 

Given Lxd = satisfied, there exists a Noether symmetry. Combined with the Lagrangian, 
this symmetry gives out 

a + 2a^ = 0, (17) 
da 

da » r a 2 dd „ 
6^-6^—^ = 0, (18) 

9/3 , . 

3a + 2a-- = 0, 19 

dcp 

3V{(f>)a + aV^((f>)f3 = 0, (20) 

which we call Noether conditions. The difference between the N-quintessence and N- 
phantom is in Eq. (jl8l) as the parameter e denotes. 

What is more, the momentum potential can be defined as below 

Pa = ^r- = 6aa, (21) 



dd eN , ■ 
— f 

,2 



a 3 0. (22) 



d<P rn z p i 

Then we can express the constant of motion which is reproduced by the Noether symmetry 

apa + fiP4> = Q = Po, (23) 

where Q is called conserved charge and is the related constant. The Noether constant 
of motion on shell gives a possibility of solving the system. More specifically, a symmetry 
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exists if at least one of the functions a or 8 is different from zero. As a byproduct, the form 
of V((f>) is determined in correspondence with such a symmetry. 

The cyclic variable can be regarded as a helpful tool of getting the exact description 
about the dynamical system. A point transformation (a, <p) — > (z, w) is effective to find the 
cyclic variable. It is 

Oz Oz 

i x z = a— + (3-j = 1, (24) 
oa d(p 

dw „dw „ . „. 

%XVJ = oi— + P—t- = 0, 25 

then the Lagrangian could be rewritten in term of the cyclic variables. After the transfor- 
mation, the cyclic variable is z, and the constant of motion can be rewritten as Q = p z . 
This will simplify our calculation effectively. A general discussion of this issue could be 
found in [TT| [T2"l ITS] . After introducing the Noether symmetry approach, we will discuss 
the solutions for the Noether conditions both in the N-quintessence and N-phantom in the 
following. 



3.1 Exact solutions for N-quintessence 

In the N-quintessence case where the sign of the kinetic terms takes the value e = 1, the 
Noether conditions are 

a + 2a^ = 0, (26) 
oa 

da „ T a 2 88 , . 

d(f> wr pl oa 
OB 

3a + 2a-- = 0, (28) 

0(p 

3aV{(f>) + aBV^cfy = 0. (29) 



When N — 1, the Noether conditions reduce to the single field case [12]. As indicated by 
Eq.(f2~Tl). the effects of the multiple scalar fields are manifested by the number of the scalar 
fields N. 

An obvious constant potential solution is 

a = 0, 8 = constant, V = constant. (30) 
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In this solution, a is the cyclic variable. And the subsequent constant of motion gives out 



Pp<t> = --^-a 3 = Q = hq. (31) 
m P i 

The discussions could be divided into two cases simply. Firstly, when /x = 0, = constant, 
this is a cosmological constant solution. Secondly, when //q ^ 0, the kinetic term <fi oc a 3 . 
The scalar fields decay fast, even faster than the corresponding vector field solution [10J. 
These two cases are trivial respectively. In the following, we will concentrate our discussions 
on another solution which is 

« = -7=> P = o 7= ' 32 
A/a za-y/a 

V = V 0(T 2 = ^ (A 2 e 2A * + B 2 e~ 2 *+ - 2AB), (33) 



where a± = Ae X(f} ± -Be A = ^JSN/Sm 2 ^ A and 5 are constants. 

We can see that the potential is a combined exponential function. Indeed, there are some 
physical origins about this kind of potential. In higher-dimensional gravitational theories 
such as superstring and Kaluza-Klein theories [14] . exponential potentials often appear 
from the curvature of internal spaces associated with the geometry of extra dimensions [15J. 
Moreover, it is known that exponential potential can arise in gaugino condensation as a non- 
perturbative effect and in the presence of supergravity corrections to global supersymmetric 
theories [16J. However, this kind of potential is picked up by Noether symmetry. 

In particular, when A = 0, the Noether conditions show 




V(<P) = V exp\-\Hr—l- ( 34 ) 



This kind of potential leads to a power-law expanding universe, with a oc t 4//3 , oj^ = —1/2. 
The quintessence with an exponential potential was widely studied in cosmology, see, for 
example, Ref. [17] . It even has a scaling solution. In the following, based on the value of 
A, we will get the exact solutions from the point of view of Noether symmetry. 

3.1.1 when A ^ and B ^ 

If we put Eqs. fj2T]) . (l22|) . ( 132]) and fp]) into Eq. (EH]), we find that the constant of motion 
is hard to obtain. Therefore, we search the cyclic variable for help. By calculating Eqs. 
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24l) and ( )25l) . we can get the following expressions for the new variables 



a 3 /V 



a 1 a + 
z = — — — , w 



6AB ' ~ 6AB ' 
where z is the cyclic variable. Correspondingly, and a could be expressed as 

, 1 . Z + W r 2 2M1/3 

0=— In , a = \9AB(z - w 2 ) 7 . 

2A 2! — w L J 

The resulting forms of potential and Lagrangian are 

4w 2 



V((/>) = V 



z 2 _ w 2 ' 



C+x = 12AB 



I -2 -2\ , 3iVVo 2 



-w 



m 



P l 



(35) 
(36) 

(37) 
(38) 



Using the Euler-Lagrangian equations, the above Lagrangian leads to the equations of 
motion for z and w, 



z — 0, w 



3NVn 



-w. 



ni 



(39) 



P i 



The solutions are 



Z — Z\t + Zq, 



I3NV 

w = wi sm{ x j — —t + w ), 



m 



(40) 
(41) 



P l 



where z , zi, wo, w\ are constants. Therefore, the exact evolution of the field and the scale 
factor could be given out as below 



Zl t + z + Wl sin(,/^t + w ) 



2A 



In 



z x t + z - [ Wl sm(J^t + w )] 

pi 



9AB((z l t + z )' 2 - wl sin 2 ' 



l3NV n 



m 



t + w )) 



P l 



1/3 



(42) 



(43) 



If z <C w, we could not get a physical value of <fi, through the scale factor seems oscillate. 
And if z ^> w, <fi is very small, but the universe will evolve as a oc t 2 ^ 3 . It is similar to the 
matter-dominated phase. 
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3.1.2 when A = and B ^ 



In this subsection, we continue to search the cyclic variables but for a different potential 
where A = while B ^ 0. By calculating Eqs. (I24j) and ( 1251) . the expressions of the new 
variables are 



7 3 / 2 



7 3 / 2 



W 



3cr + ' 3<T_|_ ' 

where z is the cyclic variable. Then <p and a can be rewritten as 

, 1 , ,B 2 ZW. /^M/a 

^=2A ln( — »- a =y ■ 

As a result, we get the potential and the Lagrangian in term of z and w 

V 



v{4>) 



Qzw^ 
-Azw N Vn 



£<t>i — — — - + 



m 2 ; 9w 2 



3 u> 2 

Apply the new Lagrangian to the Euler-Lagrangian equations, we obtain 

W n 2w 2 



2m 2 pl w' 



w 



w 



(44) 

(45) 

(46) 
(47) 

(48) 



They lead to 



w 



V Q w 2 , 3 3V w 3 2 3V w A 

t + — — — t + — — — t + w 5 



4mJ 



Am 2 pl 



4m 2 pl 



(49) 

— , (50) 

where w 2 , w 3 , w 4 are constants. Putting the above equations into Eq. (44), the evolutions 
of a and 6 are 



VgW2 ^3 _|_ 3y w 3 42 _j_ 3Vbw 4 



m v i In ( 5 2 



4m 



pi 



4m 



pi 



4m; 



pi 



W 2 £ + ^3 



(w 2 t + W 3 )( 



VoW 2 , 3 3^0^3,2 3VolU 4 



1 1/3 



4mJ 



4mJ 



4mJ 



-t + io 5 ) 



(51) 



(52) 



When z cc t, the scale factor evolves as a oc t 2 / 3 which is similar to the matter-dominated 
phase. When z ex t 3 , the scale factor is a oc t 4//3 which may accelerate the universe. This is 
an interesting solution that we need. We will discuss this solution in section fl5]) in detail. 
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However, the N-quintessence scenario could be replaced by a single field paradigm with 
the similar evolutions a oc t 4 / 3 . We just need to change the corresponding parameter in the 
single field case as 

V 0s = NV , X s = -A=. (53) 
V -/V 

The reason for this behavior is that each field experiences the 'downhill' force from its 
own potential, it feels the friction from all the scalar fields via their contribution to the 
expansion rate. 

The case B = 0, A ^ is treated exactly in the same way and the results are the same, 
except for the substitution of A for B. In summary, it must be noted that our results include 
some already known models. The exponential potential not only make the acceleration last 
a long time, but also satisfy the Noether conditions. 

As for the comparison with the observations, one field results have been derived by Ref. 
[llj. In the N-quintessence case, the range of parameter will be changed because of N. 
Considering our purpose is on the choice of the potential, we will not discuss this subject 
in detail. 



3.2 Exact solutions for N-Phantom 

For the N-phantom case where e = — 1, the Noether conditions are 

da , , 

a + 2a— = 0, 54 
da 

6^ + iv4f = 0, (55) 
o<p oa 

d/3 , . 

3a + 2a-- = 0, 56 

ocp 

3aV(<P) + af3V; ht (<f ) ) = 0. (57) 

Compared to the N-quintessence case, the difference arises in Eq. (|55p by the sign of the 
kinetic terms. 

Obviously, the simplest solution is 

a = 0, (3 = constant, V = constant. (58) 

This constant potential solution is similar to the corresponding solution in the N-quintessence 
case, we don't discuss this fast decaying case. 
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However, another interesting solution is 



2Ccos(±,/^^M 

V2 V 2 m pl' /r n \ 

a = -= , (59) 

>a 



= / 8 ^ , (60) 

aJa 



2 1 j3N JN V Q „ /3N 



V(<P) = Voting -T+) = f^ 1 - COS (V (61) 

where C is a constant. When N — 1, there are some differences between the results in Ref. 
[18J and ours. The form of the potential could be called PNGB (Pseudo Nambu-Goldstone 
Bosons) potential resulting from explicit breaking of a shift symmetry [19]. 

To find the exact evolution of the universe, as the calculations in the N-quintessence case, 
we need the help of the cyclic variables. According to Eqs. ( J24|) and (1251) . a transformation 
could be done from (a, 0) to (z,w), 

= arctan -, a = (3C) 2/3 (z 2 + w 2 ) 1/3 } (62) 

z 



then we can rewrite the potential and the Lagrangian as 

V w 2 
z 2 + w 



(63) 



£*i = 9C 2 



4.. 2 , 2 , Nw 
-(z 2 + w 2 )+V — 
3 m£ 



(64) 



The Lagrangian leads to the equations of motion for the new variables 

3NV w 
4 m 



z = 0, w = ^ — • (65) 



pi 

The solutions are 

z = z 3 t + z 2 , (66) 

w = io 6 exp(W~ — 2-t), (67) 
V 4 "V 

where Z2, £3 and wq are constant. However, by using the cyclic variable z, we get the 
evolutions of the field and the scale factor, 



»eexp(Jft) 

arctan — , (68) 

z 3 t + z 2 
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(3C) 



2/3 



(z 3 t 



Z2) 2 + w^expi 



<3NVn 



1 1/3 



ni 



t) 



P l 



(69) 



If z 3> w, the values of fields are nearly zero, a oc it is the matter-dominated solution. 
When z <C w, the universe evolves as a oc exp( 
will discuss this solution in section (151) in detail 



3A'\„ 



T^t), this is the de-Sitter solution. We 



As for the comparison with the observations, one field results have been derived by 
Ref. [IB]. In the N-phantom case, the range of parameter will be changed because of N. 
Considering our purpose is on the choice of the potential, we will not discuss this subject 
in detail. 



4 Noether Symmetry in N-Quintom Case 

The quintom scenario is proposed to fit the observable data [5]. N-quintessence and N- 
phantom could not cross = —1 as we see. However, N-quintom has an attractive feature 
that it may cross = — 1 which is a possibility implied by the data. After adding the 
Noether symmetry, this property should be rechecked. Though in the "cosmic triad" vector 
field case, Noether symmetry provides an interesting constraint on the potentials [10] for the 
quintom case with crossing —1. However, it is worthy of trying the Noether symmetry 
approach in the N-quintom case. According to the action (j5]), the point-like Lagrangian is 



Co 



3aa 2 




q - - a 3 V r , 




C^Vp + PmO- 



(70) 



Based on the above point-like Lagrangian, the total energy and the Euler-Lagrangian 
equation will give out the Friedmann equation, the Raychaudhuri equation and the equa- 
tions of motion 



H 2 



H 



1 



m 



pi 



N q (f + V q ) + N p (-^ + V p )+p r 



b q + 3H<j> q + VI 

b v + 3H<p p - v; 



+ P<t> + Pr, 

0, 
0. 



2 q 2 p 



m 



pi 



(71) 

(72) 

(73) 
(74) 



where the primes mean VL = dV q /d(j) q and VL = dV p /d<p p . The energy density and the 
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pressure which could be derived from the action (jSJ) are 



•2 • 2 

= N q (^ + K 9 ) + iV p (-^- + V p ), (75) 

; 2 -2 
= N q (^- - V q ) - N p (^ + V p ). (76) 

(77) 

So the EoS parameter is 

— 2 j- 2 • (78 J 

N q (*t + V q ) + N p (-%- + V p ) 

Now, we should choose a new configuration space Q = (a, <f) q , <f> p ) with the correspond- 
ing tangent space TQ = (a, (f> q , (f> p , a, (f> q , <f> p ). And the vector generator which induce the 
Noether symmetry is changed to 

9 s d d ^ d ~ d d . . 

X = «— + + 7 — + «— + /3- r + 7 - r , 79 

da d(f) q d(j) p da d(j) q d<p p 

where 5, B and 7 are generic functions of the variables a, cj) q and <f> p . The Noether symmetry 
requires the Lie derivative of the Lagrangian vanishes which means -^£2 — 0. Following 
Ref. [TTJ [TSJ [32], the Noether conditions can be obtained 

da , , 

5 + 2a— = 0, 80 
da 

6^-^ = 0, (81) 
d(p q m z pl da 

Q da_ + Npa^dj = Q ^ 
<90 p m 2 pl da 

88 , . 

3a + 2a-— = 0, 83 
d(p q 

35 + 2a^- = 0, (84) 
d(p p 

S(V q + V v )a + aV^~B + aV^y = 0. (85) 
There is an obvious solution that is 

5 = 0, = D, j = D, (86) 
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where D and D are integral constants. The symmetry exists, if and only if at least one of 
the parameter a, (3, 7 is not zero. Based on the Noether conditions, we find a condition 
relating the potential forms of the quintessence-like fields and the phantom-like fields, that 
is 

DV; q = -DV; p . (87) 
And the constant of motion corresponding to this solution is 

- DN q a 3 <p q + DN p a 3 <p p = Q = /i . (88) 

In the following discussion, based on the value of D, D and /i , we try to discuss the 
solutions, especially for the value of EoS parameter. 

4.1 When D ^ and D = 

If D 7^ and D = 0, we can get VI = 0, the quintessence-like matter has a constant 
potential. And from the constant of motion, we can get <p q oc /ipa -6 . However, based on 
the value of /x , we divide the situation into two cases to discuss. 

Case a), when fiQ 7^ 0, the kinetic terms of the quintessence decay fast, while their 
potentials are constant, and no constraint on the phantom type of matter, which leads to 
w < —1 at last. 

Case b), when /xq = 0, the quintessence scalar filed is a constant. This case is similar to 
a phantom model with cosmo logical constant. The interesting thing is that we could not 
give any constraint on the phantom-like matters. 

The case D = and D^O could be treated exactly in the same way. And the results 
are the same, except for the non-constrained field is changed to the quintessence-like type. 

4.2 When D ^ and D ^ 
4.2.1 The fj, ^ case 

In this case, the conserved charge is not zero. From Eq. (188]) . we get that <p q = D<p p /D oc 
a -3 . It means that the kinetic terms of the scalar field decay fast. The equations of motion 
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leads to VL = Vi = 0, i.e., the potentials are constant. However, the EoS parameter 
evolves to = — 1 until the kinetic terms of the scalar fields vanish. 

4.2.2 The /i = case 

In this case, the conserved charge vanishes, so DN q <p q = DN p <fi p , combined with Eq. ( 1871) 
and the equations of motion, N q = N p is obtained. We put these results into Eq. (17H|) . and 
get 



q --v q -v p 



U(f) = 1— » *. (89) 



2 h P 

: 2 K2 



If D/D < 1, (f) p < (f) q) > — 1. The physical meaning is that if the quintessence type 
fields slowly vary compared with the phantom type fields, the quintessence will take the 
dominating role, and make > —1. And we can discuss the D/D > 1 case in the same 
way, where the phantom type fields will take the dominating role and < —1. However, 
this solution is new. And if it cross = —1, the ratio D/D should be variable. However, 
as Noether symmetry approach required, D/D is constant. It means in N-quintom case, 
after adding Noether symmetry, we could not make this scenario cross = —1. 

In a short summary, even the Noether symmetry does not give an explicit potential 
in N-quintom case, it gives a constraint on the forms of the scalar field potentials. If we 
try to connect this model to the observations such as SNIa data, we must choose a proper 
potential. Unfortunately, the observations will give constraints to the potential parameter 
not the parameter related to Noether symmetry which we are interested here. And this 
symmetry restricts the EoS parameter of crossing = —1. 



5 From Deceleration to Acceleration 

Based on the exact potential forms given by Noether symmetry in N-quintessence and 
N-phantom models, the evolution of our universe could be analyzed. Firstly, two new 
variables y = <f)/m p i, u = ln(a/ao) are needed. Then we can define the fractional energy 
density of dust matter as fl m = p m /?>H 2 m'^ l = Q m0 (H /H) 2 exp(— 3m), and the fractional 
energy density of scalar fields fi^ = p^/SH 2 ™^ which depends on the exact potential form. 

In N-quintessence model, we discuss the possible accelerating solution which is presented 
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Figure 1: The evolutions of fractional energy densities and VL m in N-quintessence model. 



in Eq. ( 133]) with ^ = 0,5^0 and A = \/3/2. can be written down as 



,'2 



fi^ = — + IV exp(-Ay), 
6 

where fiy = V B 2 /3H 2 m 2 pl = Q V0 {H /H) 2 . Then, we can simplify Eq.® and (HU) as 

# 2 fi m0 exp(-3«) + fiy exp(-Ay) 



(90) 



y 



3Afiy exp(— Ay) 



1 - y' 2 /6 
3 



f2 m + 3f2y exp(— Ay) 



(91) 
(92) 



where a prime denotes the derivative with respect to u. Following the numerical calculation 
method used in Ref . [TH [20] , the evolution of the fractional energy densities can be plotted. 
We choose Q m = 1 in the matter dominated epoch around a/a^ ~ 1/12 or u = —2.5 as 
initial condition. FigJTJ shows today's fractional density density f2 m0 is nearly 0.6 which is 
contradictable with the widest observational results Q m0 = 0.3 ±0.1 
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Figure 2: The evolutions of fractional energy densities and Vt m in N-phantom model. 



Furthermore, setting N = 1 and using the potential in Eq. (l6Tj) . £1$ in N-phantom case 
reads 

a, — • SMI -ci^h/— //)). (93) 



— + n v (l - cos{\l —y)), 



where fly — V /3H 2 m 2 l = Q VQ (H /H) 2 . Then, the evolutions of scale factor and scalar 
field in N-phantom case are 



H 2 fl m0 exp(-3u) + Q V0 (1- cos(y^f y)) 

( #o } = l - y' 2 /Q : 



V = --Afiysm(W — y) - 



-il m + — — (1 - cos( 



3N 



y)) 



y 



(94) 
(95) 



We can also start from the matter dominated epoch around a/ao ~ 1/12 oi u = —2.5, and 
give out the evolutions of the fractional energy densities. Figj2] shows today's fractional 
density of dust matter Q m Q is nearly 0.23 which is consistent with the observational results 
Q m0 = 0.3 ±0.1. 
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Figure 3: The evolutions of the minus of the deceleration factor — q in N-quintessence and 
N-phantom cases. 



Specifically speaking, we can write down the acceleration (the minus of the deceleration 
factor) 

^ = - q = n^- l - y * - Ia m , (96) 

and plot its evolutions in N-quintessence and N-phantom cases. FigfJ shows N-quintessence 
with exponential potential chosen by Noether symmetry cannot make our universe acceler- 
ate (— q < 0), while N-phantom with the cosine potential can make our universe accelerate 
(— q > 0). These results coincide with the evolutions of the fractional energy densities in 
FigJT]and|2j In conclusion, today's acceleration heavily depends on the choice of potentials. 
This is the reason why we don't discuss the evolution of our universe in the N-quintom case 
where the Noether symmetry doesn't give out the exact form of potentials. 
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6 Conclusion 



There is no immediate physical justification for the choice of V(4>) in multiple scalar fields. 
In this Letter, to choose proper potentials for multiple scalar fields scenario, and to be 
consistent with the observations which indicates the EoS parameter in the range of —1.21 < 
u<f, < —0.89, we have studied the N-quintessence, N-phantom, N-quintom scalar fields 
models by the Noether symmetry approach. The existence of Noether symmetry implies 
that with respect to the infinitesimal generator of the desired symmetry, the Lie derivative 
of the related Lagrangian vanishes. As we have considered a flat FRW metric, the phase 
space in the N-quintessence and N-phantom was then constructed by taking the scale factor 
a and the scalar field as independent dynamical variables. In the N-quintom case, we 
have to expand the configuration space to Q = (a, <f> q , 4> p ). 

Specifically speaking, on the one hand, the Noether conditions depend on the cosmo- 
logical dynamics which is determined by the potentials. On the other hand, the main con- 
sequence by adding the Noether symmetry is that we have selected the class of potentials 
and indicated the most reasonable, specific ones directly from the physical interpretation. 
In the N-quintessence case, we find the exponential potentials from the Noether conditions 
which could be derived from several theoretic models. In the N-phantom case, the suitable 
potential required by the Noether conditions is -y(l — C0S (\J^^^~ l )) which is related to 
pseudo Nambu-Goldstone boson. The case of the N-quintom is very interesting. Although 
it does not give an explicit potential, it gives a constraint on the forms of the scalar field 
potentials. 
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